Singular Structure and Enhanced Priedel Oscillations in the Two-Dimensional 

Electron Gas 
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We calculate the leading order corrections (in rg) to the static polarization n*(q, 0, ), with dynam- 
ically screened interactions, for the two-dimensional electron gas. The corresponding diagrams all 
exhibit singular logarithmic behavior in their derivatives at g = 2kF and provide significant enhance- 
ment to the proper polarization particularly at low densities. At a density of = 3, the contribution 
from the leading order fluctuational diagrams exceeds both the zeroth order (Lindhard) response 
and the self-energy and exchange contributions. We comment on the importance of these diagrams 
in two- dimensions and make comparisons to an equivalent three-dimensional electron gas; we also 
consider the impact these finding have on 11* (q, 0) computed to all orders in perturbation theory. 



Two-dimensional electron systems with standard 
singular Coulombic interactions have attracted a great 
deal of attention in recent years. Of particular interest 
have been the driving mechanisms behind the high tem- 
perature superconductors, the metal-insulator transition 
exhibited by the two-dimensional localized electrons in 
semiconductors, and the novel charge and spin-density- 
wave ordering found in layered compounds. An impor- 
tant question is the degree to which we might attribute 
this novel behavior specifically to the two-dimensionality 
and the singular character of the Coulomb interaction? 
In this letter, we attempt to provide an answer by com- 
puting the leading order corrections to the proper static 
polarization 11* beyond the random-phase approximation 
(RPA) through inclusion of the diagrams presented in 
Fig. shown there as the self-energy, exchange, and 
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FIG. ^ (a) Leading order corrections to the proper polar- 
ization n*. The interaction lines represent RPA screened 
Coulomb interactions, (b) The three-point diagram A'^^'and 
the effective scattering vertex r(q) . 



fluctuation diagrams, respectively. For the three- 
dimensional electron gas, these diagrams have been well 
studied, and in particular have been computed and ana- 
lyzed in detail (see 0- and references therein). For a 
strictly two-dimensional electron gas with 27re'^ /q inter- 
action, we report here that the response function beyond 
the RPA exhibits singular structure that is far more pro- 
nounced than found from the zeroth order Lindhard re- 
sponse function; these leading order corrections play an 
important role even at high-densities. 

It is not surprising that all three leading order di- 
agrams exhibit similar singular structure since they can 
be reduced to expressions containing the three-point di- 
agram shown A'^'^Hn Fig. [| (multiplied by Coulombic in- 
teraction lines). This singular structure should extend 
to all orders in perturbation theory since all higher or- 
der corrections to the polarization also reduce to a sum of 
three-point diagrams A*^'^^ (again multiplied by Coulombic 
interaction lines) . More importantly, the enhancement of 
the polarization from all of three of the leading order di- 
agrams is already significant at high densities, but as the 
density is reduced the contribution from the fluctuation 
diagrams rapidly exceeds that of the self-energy and ex- 
change contributions. This is in marked contrast to the 
three-dimensional electron gas where fluctuational dia- 
grams provide a small enhancement to the correspond- 
ing zeroth order Lindhard response when compared to 
contributions solely from the exchange and self-energy 
diagrams (for densities in the metallic range |^). The 
increased role that fluctuations are expected to play in 
two dimensions is evident. 

Begin with the static zeroth-order Lindhard re- 
sponse function 11'^, appropriate to the random-phase 
approximation. The well-known singular behavior of 
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n" leads to possible instabilities arising from electron- 
electron interactions. In three dimensions, this singular- 
ity also leads to a modulation in the long range effective 
electron-electron interaction which in three dimensions 
takes the form V{r) ~ cos{2kFr)/r^ . Based on this Kohn 
and Luttinger Q| were the first to consider the possibility 
of ground-state electron pairing in the presence of purely 
repulsive forces where they found attractions in the ef- 
fective vertex T[q) of Fig. ^for large angular momentum 
quantum numbers. Details of the shape of the Fermi 
surface (for example, nesting) might well enhance such 
possibilities. 

In two dimensions, and at the level of the Lindhard 
response function, these effects are not especially promi- 
nent. The function n"(g) is then a constant for q < 2fci?, 
and exhibits a square root singularity in the derivative 
for momentum q > 2kp Because of this arguments 
for possible charge/spin instabilities based on are not 
compelling. While this square root singularity does lead 
to a long range component in the effective interaction 
(y(r) ^ sin{2kpr)/r'^) for two dimensions, the effective 
interaction vertex T{q) of Fig. |] computed at this low or- 
der remains independent of momentum at the Fermi sur- 
face. Only an s-wave repulsive scattering channel opens 
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at this order and higher-order effects must be considered 
in order to find attractive scattering channels Q. But 
beyond the random phase approximation (RPA) the sit- 
uation is quite different and we report here that the first 
order corrections to 11* already give singularities (in part 
logarithmic) in the derivative on both sides of 2kp, and 
this singular structure is further enhanced according to 
the number of three-point diagrams A^^-'and interaction 
lines appearing in a polarization diagram. As an exam- 
ple, the fluctuation diagrams 11// of Fig. [l] contains two 
three-point diagrams and two dynamically screened in- 
teraction lines, and shows greater singular behavior at 
the Fermi surface in two-dimensions than the self-energy 
and exchange diagrams which depend on a single three- 
point diagram and single dynamically screened interac- 
tion line (see Eqs. |^ and ||). Consequently the enhance- 
ment near 2k p in 11* from these diagrams leads to a 
greater attraction in real space for the screened interac- 
tion V* = f(j/(l — VqQ*) when compared with the RPA 
effective interaction V^^^ = Vq/ {1 — VqUP). 

As noted in references and [Q, the self energy 
Ilse , exchange YL^x , and fluctuation 11// diagrams de- 
picted in Fig. |l| can all be written in terms of the three- 
point function; thus 
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Here p — {iujp, p) is the energy-momentum variable, 
vrpa is the screened Coulombic interaction in the ran- 
dom phase approximation (using Vq = 27re^/g), and Tvq 
stands for the two-dimensional trace over all momentum 
and energy variables, i.e. / dujq/{2TT) J Cpq/ {2-k)'^ . The 
requisite three-point function A(^-'(g,p) , is given by 

A^^H?,?) = -2Tr^^Go{k)Go{k + q)G^{k + q + p), (4) 

where exact analytical expressions for two-dimensional 
forms are given by Neumayr and Metzner We com- 
pute the diagrams of Fig. ^ in imaginary frequency space 
to circumvent nonintegrable divergences near 2kp and 
also when integrating through plasmon peaks in the 
screened interaction VRPA^q)- To obtain the static po- 
larizability we simply perform the requisite analytic con- 
tinuation n(p, iujp) — > n(p, Lu + i5) = n(p, o + iS). 

The numerical procedure for computing Ilse , ^ex , 
and 11// is lengthy but straightforward. The results for 
a density of Ts = 2 (not a particularly high density) 
are plotted in Fig. ^ and can be summarized as follows: 
First, Hex exceeds 11° for momenta up to and slightly 



above 2kp, but Ilse and H^x have opposite signs and 
their sum largely cancels. However, H// at this density, 
already gives a slightly larger contribution than H^e and 
Hex combined. Second, all three diagrams are finite but 
discontinuous in slope at 2kp; their derivatives diverge 
logarithmically at all densities. We have determined 
through a numerical fit that the logarithmic structure 
is of the form y^\q — 2kp \ log |g — 2kp\ on either side of 
q = 2k p. Maldague has numerically computed the 2D 
self-energy and exchange contributions in the absence of 
screening and found similar logarithmic divergence in the 
slope for q > 2kp by invoking an electrostatic analogy 
that treats the integrals over {Use + ^ex) as equivalent 
charge distributions of radius kp concentrated near the 
origin times the 2D bare interaction. The divergence in 
the slope is stronger above 2kp in overall magnitude than 
below 2kp due to the analytic form of A^^ Which contains 
non-trivial expressions involving square root singularities 
(similar to that found in H" ). These expressions con- 
tribute to the analytic form of A^^^ above, but not below 
2kp. Overall the divergent structure is greatest for H// as 
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stated earlier. The contribution from all three diagrams 
is significant already at = 2. The role of fluctuations 



11// can be thought of as representing quantum fluctua- 
tions in the screening clouds surrounding two interacting 
electrons, and we expect such effects to play a greater 
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FIG. The solid curve is the static 11/; , the dashed curve 
shows {Use + Hex), and the dashed-dotted curve is the sum of 
all three of the first order correction all at = 2. The inset 
contains the fluctuation diagram at a lower density of Vs = 4. 
The dotted curve is the 2D Lindhard response. 

is expected to progressively dominate with decreasing 
density, and this is apparent at rj = 3 where the fluc- 
tuation diagrams begin to exceed 11" . By = 8 (a 
typical density for high-temperature superconductors or 
the two-dimensional electron(hole) gases in MOSFETS) 
it is a factor of four greater than 11° where at this den- 
sity the sum of the self-energy and exchange diagrams 
only just begins to exceed the zeroth order Lindhard re- 
sponse. Overall the effect of dynamical screening as 
increases is to reduce the total enhancement of the polar- 
ization diagrams as well as the degree of divergence in the 
slope (see Q for a comparison with the self-energy and 
exchange diagrams computed with the bare interaction) . 

The leading order corrections to the polarization be- 
yond the RPA are not as significant in three-dimensions. 
In the absence of screening the enhancement (Ilse -|- 
Ilex)/!!" is O.ITts at p = compared with OAbvs in 
2D , and then proceeds to sharply falls off to half its 
value at 2A:f The singular structure found in these 
diagrams resembles that of the 3D Lindhard response, 
but with dynamical screening the overall enhancement 
of these diagrams is greatly reduced as shown in Fig. ^ 
(the data is from reference |^). This data is represen- 
tative of the entire metallic density range and confirms 
the expectation that self-energy and exchange effects are 
more pronounced in 2D. More striking is the significant 
enhancement of n*^ in 2D from the fluctuation diagrams 
as seen by comparing the plots of Fig. ^ with the three- 
dimensional data of Fig. ^ Physically, the diagrams of 
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FIG. |. The Static -Use +nex , and Ufi in 3D at Ts = 2. The 
inset shows the 3D Lindhard response. 



role in 2D |1C|J11[]. Enhanced fluctuations in two- 
dimensions are responsible for the destruction of super- 
conducting phase coherence and long range crystalline 
order for T > 0, and thus it is not surprising to see this 
theme played out here as well 
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As mentioned earlier, the general structure of these 
leading order diagrams is an integral over three-point 
functions A augmented by Coulombic interaction lines. 
For higher orders, an explicit reduction formula derived 
by Neumayr and Metzner [^j holds that the general N- 
loop diagram can be expressed in terms of 3-loop dia- 
grams over an appropriate energy/momentum transfer 
factor (i.e. in Il^x that factor is (p-q) in Eq.(^) and thus 
the structure of all higher order diagrams will be of sim- 
ilar form (for an explicit formula see reference Q). We 
therefore expect the singular logarithmic features found 
in these lower order diagrams to be present in all higher 
order corrections. Cancellations between some of the 
higher order diagrams may well reduce the overall sin- 
gular behavior of the total proper polarization H* . 

As further higher order corrections are considered 
there will continue to be significant cancellation between 
vertex and self-energy diagrams in such a way that their 
sum makes a smaller net contribution, and any reason- 
able higher order calculation to H* should include the 
right mix of vertex and self-energy diagrams (as is the case 
for the conserving approximation of Baym and Kadanoff, 
for example) . But the contribution from diagrams of 
the fluctuation type differs in two respects. First, the 
internal Coulomb lines are screened and already repre- 
sent the sum of an entire class of diagrams in a Dyson 
series sense Gffl; second, higher order contributions of 
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diagrams from the fluctuational type do not cancel but 
actually sum to contribute to quantities such as the pair 
correlation function at zero separation, g{0) We 
generally expect, therefore, that contributions from all 




FIG. g| The screened effective interaction in real space show- 
ing an enhancement in the amplitude of the Friedel Oscilla- 
tions due to the singular structure in 11* at rs = 2 (the dashed 
line is the RPA screened interaction). The inset shows the 
three dimensional case at an equivalent density. 

such diagrams (e.g. the ladder series) will significantly 
contribute to the total proper polarization H* leading 
to an overall enhancement and to singular logarithmic 
structure at 2k p. 

Finally, we return to the important Kohn-Luttinger 
question for the two dimensional case, but in the context 
of the leading order corrections of Fig. ^ The effective in- 
teraction vertex T{q) in the case of the electron gas must 
include a sum of all polarization diagrams due to the sin- 
gularity at small momentum transfer q = p — p'. If we 
include the leading order contributions to T{q) we obtain 
T{q)=vJ{l-VgU*{q))with U* = U° + Use + Il,, + Ufi. 
The logarithmic singularity in H* now gives rise to a an 
attraction in T; for even angular momentum / starting at 
I = 2 (d — wave), behavior not possible if only Xl" entered 
into the argument. These singularities also significantly 
contribute to the modulation in real space of the equiv- 
alent effective interaction (i.e. V{q) = Vq/{1 — Vqll*)) 
now plotted in Fig. ^ for rg = 2. The principal depth 
(the first minimum) of the effective interaction is strik- 
ing, and the effect is further enhanced as the fluctuation 
diagrams exceed for densities as low as rg = 3. The 
inset in Fig. ^ shows the three-dimensional case for com- 
parison where the singularity at 2k p from the leading 
order corrections leads to only a small enhancement of 
the Friedel oscillations. 

To summarize, the low order diagrams, reducible to 
similar integral forms, all exhibit singular logarithmic be- 



havior at 2kp and provide significant enhancement over 
the zeroth order Lindhard response. Real systems ei- 
ther have either a finite transverse extent or exhibit a 
multi-layered structure. To account for this, we have 
substituted approximate but more realistic forms for the 
potential in the diagrams of Fig. |^ more appropriate 
to 2D electrons in metal-oxide-semiconductors, in thin 
metallic films, or in layered superlattices. Similar re- 
sults are found for systems where realistic thicknesses 
and superlattice spacings are used except that the over- 
all enhancement from these diagrams compared to H'^ is 
reduced depending on the thickness (or spacing between 
the planes), dielectric constant, and the effective mass of 
the system. Our results extend to other two-dimensional 
systems and apply to the 2D electron or hole gas systems 
found in metal-oxide-semiconductors. The debate over 
the nature of the metallic state in the recently observed 
metal-insulator transition of 2D electron or hole semi- 
conductor systems still ensues with explanations ranging 
from superconductivity to more exotic non-Fermi liquid 
behavior. The present results may bear on (and lend sup- 
port to) explanations of the former type. More impor- 
tantly, we note that accurate data at 2kp is of consider- 
able interest, particularly from Monte Carlo simulations. 
Large error bars generally accompany polarization data 
from Monte Carlo simulations in the intermediate regime 
around 2k p and make difficult an accurate determination 
of quantities such as local field factor corrections and cor- 
responding effective interactions that could lay to rest 
debates concerning the likelihood of superconductivity 
or other electronic instabilities in these systems. 
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